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We point out a new and largely model-independent constraint on the dark matter scattering cross 
section with nucleons, applying when this quantity is larger than for typical weakly interacting 
dark matter candidates. When the dark matter capture rate in Earth is efficient, the rate of 
energy deposition by dark matter self-annihilation products would grossly exceed the measured 
heat flow of Earth. This improves the spin-independent cross section constraints by many orders of 
magnitude, and closes the window between astrophysical constraints (at very large cross sections) 
and underground detector constraints (at small cross sections). In the applicable mass range, from 
~ 1 to ~ lO'^'' GeV, the scattering cross section of dark matter with nucleons is then bounded from 
above by the latter constraints, and hence must be truly weak, as usually assumed. 

PACS numbers: 95.35.-|-d, 95.30.Cq, 91.35.Dc 



I. INTRODUCTION 

There is a large body of evidence for the existence of 
dark matter, but its basic properties - especially its mass 
and scattering cross section with nucleons - remain un- 
known. Assuming dark matter is a thermal relic of the 
early universe, weakly interacting massive particles are 
prime candidates, suggested by constraints on the dark 
matter mass and self-annihilation cross section from the 
present average mass density m. However, as this re- 
mains unproven, it is important to systematically test 
the properties of dark matter particles using only late- 
universe constraints. In 1990, Starkman, Gould, Es- 
mailzadeh, and Dimopoulos examined the possibility 
of strongly interacting dark matter, noting that it indeed 
had not been ruled out. Many authors since have ex- 
plored further constraints and candidates. In this litera- 
ture, "strongly interacting" denotes cross sections signif- 
icantly larger than those of the weak interactions; it does 
not necessarily mean via the usual strong interactions be- 
tween hadrons. We generally consider the constraints in 
the plane of dark matter mass and spin-independent 
scattering cross section with nucleons a^N- 

Figure [T] summarizes astrophysical, high-altitude bal- 
loon/rocket/satellite detector, and underground detector 
constraints in the a^N-m^ plane. Astrophysical limits 
such as the stability of the Milky Way disk constrain 
very large cross sections 0, Q • Accompanying and com- 
parable limits include those from cosmic rays and the 
cosmic microwave background 0, Small cross sec- 
tions are probed by CDMS and other underground de- 
tectors d, 0, H, @] ■ A dark matter (DM) particle can be 
directly detected if a^N is strong enough to cause a nu- 
clear recoil in the detector, but only if it is weak enough 
to allow the DM to pass through Earth to the detector. 
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FIG. 1: Excluded regions in the a-^N-m^ plane, not yet in- 
cluding the results of this paper. From top to bottorn, these 
come from astrophysical constraints (dark-shaded) 0, '3,0,1^, 
re-analyses of high-altitude detectors (medium-shaded) 2, 1^, 
[TTI . [l^ , and underground direct dark matter detectors (light- 
shaded) 0, 0, B 01 • The dark matter number density scales as 
and the scattering rates as a^N/m^\ for a fixed scat- 
tering rate, the required cross section then scales as m^. We 
will develop a constraint from Earth heating by dark matter 
annihilation to more definitively exclude the window between 
the astrophysical and underground constraints. 



In between the astrophysical and underground limits 
is the window in which a-^^N can be relatively large [2|. 
High-altitude detectors in and above the atmosphere 
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have been used to exclude moderate-to-strong values of 
the cross section in this region 0, [13, [Uj ■ However, 
there are still large gaps not excluded. There also is 
some doubt associated with these exclusions, as some of 
the experiments were not specifically designed to look for 
DM, nor were they always analyzed for this purpose by 
people associated with the projects. In fact, the exclu- 
sion from the X-ray Quantum Calorimetry experiment 
was recently reanalyzed [l3| and it changed substantially 
from earlier estimates [ll|. If this intermediate region 
can be closed, then underground detectors would set the 
upper limit on a-^jq- That would mean that these detec- 
tors are generally looking in the right cross section range 
and that DM-nucleon scattering interactions are indeed 
totally irrelevant in astrophysics. 

We investigate cross sections between the astrophysi- 
cal and underground limits, and show that u-^n is large 
enough for Earth to efficiently capture DM. Incoming 
DM will scatter off nucleons, lose energy, and become 
gravitationally captured once below Earth's escape ve- 
locity (Section 4). If this capture is maximally efhcient, 
the rate is 2 x 10^^ (GeV/m^) s^-^. The gravitationally- 
captured DM will drift to the bottom of the potential 
well. Earth's core. Self-annihilation results if the DM is 
its own antiparticle, and we assume Standard Model final 
state particles so that these products will deposit nearly 
all their energy in the core. 

Inside a region in the a^N-m^ plane that will be de- 
fined, too much heat would be produced relative to the 
actual measured value of Earth's heat flow. The max- 
imal heating rate obtained via macroscopic considera- 
tions is ~ 3330 TeraWatts (TW), and follows the max- 
imal capture rate, assuming that Earth is opaque with 
a geometric cross section. Note that the flux of DM 
scales as l/m^, while the heat energy from annihila- 
tions scales as m^, yielding a heat flow that is inde- 
pendent of DM mass. The efficient capture we con- 
sider leads to a very similar heating rate, though it 
is based on a realistic calculation of microscopic DM- 
nucleon scattering, as discussed below. DM interactions 
with Earth have been previously studied in great detail, 
e.g., Refs 0, M, [3, \m, M, [O, M, Hi, [13, [III, but those 
investigations generally considered only weak cross sec- 
tions for which capture is inefficient. 

In our analysis, the cr^jv exclusion region arises from 
the captured DM's self-annihilation energy exceeding 
Earth's internal heat flow. This region is limited below 
by the efficient capture of DM (Section 4), and above by 
being weak enough to allow sufficient time for the DM 
to drift to the core (Section 5). These two limits define 
the region in which DM heating occurs. Why is it im- 
portant? Earth's received solar energy is large, about 
170,000 TW but it is aU reflected or re-radiated. 
The internal heat flow is much less, about 44 TW (Sec- 
tion 3) [13|. Inside this bounded region for a-^N, DM 
heating would exceed the measured rate by about two 
orders of magnitude, and therefore is not allowed. We 
will show that this appears to close the window noted 



TABLE L Relevant heat flow values. The top entries are 
measured, while the lower entries are the calculated potential 
effects of dark matter. 



Heat Source 


Heating Rate 


Solar (received and returned) 
Internal (measured) 


170,000 TW 
44.2 ± 1 TW 


DM annihilation (opaque Earth) 
DM annihilation (our assumptions) 
DM kinetic heating 


3330 TW 
3260 TW 
~ 3000 X lO"*^ TW 



above in Fig. [1] up to about ~ 10^° GeV. In order 
to be certain of this, however, we call for new analyses 
of the aforementioned constraints, especially the exact 
region excluded by CDMS and other underground detec- 
tors. Our emphasis is not on further debate of the details 
of speciflc open gaps, but rather on providing a new and 
independent constraint. In Table I, we summarize the 
heat values relevant to this paper. While the origin of 
Earth's heat flow is not completely understood, we em- 
phasize that we are not trying to account for any portion 
of it with heating from DM. 

There has been some previous work on the hea ting 
of planets by DM annihilation [H, [H, [H, [H, [13, [M 
[29, 30]. These papers have mostly focused on the Jo- 
vian planets, for which the internal heat flow values 
are deduced from their infrared radiation (sij . In some 
cases [il,[li,[li,[li,[13, DM annihilation was invoked to 
explain the anomalously large heat flow values of Jupiter 
and Saturn, while in other cases [H, [13, [sOj , the low heat 
flow value of Uranus was used to constrain DM annihi- 
lation. An additional reason for the focus on these large 
planets is that they will be able to stop DM particles of 
smaller cross section than Earth can (Ref. [23] consid- 
ered Earth, but invoked an extreme DM clumping factor 
to overcome the weakly interacting cross section). How- 
ever, as we argue below in Section 6, the more relevant 
criterion is how significant of an excess heat flow could 
be produced by DM annihilation, and this is much more 
favorable for Earth. (If this criterion is met, then the 
ranges of excluded cross sections will simply shift for dif- 
ferent planets.) Furthermore, the detailed knowledge of 
Earth's properties gives much more robust results. In this 
paper, we are presenting the first detailed and systematic 
study of the broad exclusion region in the a^jq-m^ plane 
that is based on not overheating Earth. 

Our constraints depend on DM being its own antipar- 
ticle, so that annihilation may occur (or, if it is not, that 
the DM-antiDM asymmetry not be too large). This is 
a mild and common assumption. The heating due to 
kinetic energy transfer is negligible. Since the DM ve- 
locity is ~ 10~'^c, kinetic heating is ~ 10~^ that from 
annihilation, and would provide no constraint (Section 
4). The model-independent nature of our annihilation 
constraints arises from the nearly complete insensitivity 
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to which Standard Model particles are produced in the 
DM annihilations, and at what energies. All final states 
except neutrinos will deposit all of their energy in Earth's 
core. (Above about 100 TeV, neutrinos will, too.) Since 
the possible heating rate (> 3000 TW) is so large com- 
pared to the measured rate (~ 40 TW), in effect we only 
require that not more than ~99% of the energy goes into 
low-energy neutrinos, which is an extremely modest as- 
sumption. 

Some of the annihilation products will likely be neu- 
trinos, and these may initiate signals in neutrino detec- 
tors, e.g., as upward-going muons 0, [H, [H, HI, [H, [3^, 
113, [ssl. l39l|. While the derived cross section limits can 
be constraining, they strongly depend on the branching 
ratio to neutrinos and the neutrino energies. Compre- 
hensive constraints based on neutrino fluxes for the full 
range of DM masses appear to be unavailable; most pa- 
pers have concentrated on the 1-1000 GeV range, and a 
few have considered masses above 10* GeV. We note that 
the constraints for DM masses above about 10^" GeV 
may require annihilation cross sections above the unitar- 
ity bound, as discussed below. As this paper is meant to 
be a model-independent, direct approach to DM proper- 
ties based on the DM density alone, we do not include 
these neutrino constraints. 

We review the current DM constraints in Section 2, 
review Earth's heat flow in Section 3, calculate the DM 
capture, annihilation, and heating rates in Sections 4 and 
5, and close with discussions and conclusions in Section 
6. 



II. REVIEW OF PRIOR CONSTRAINTS 

Figure [1] shows the current constraints in the ct^at- 
plane. As we will show, the derived exclusion region 
found by the requirement of not overheating Earth using 
DM annihilation lies in the uncertain intermediate area 
between the astrophysical and underground constraints. 

A. Indirect Astrophysical Constraints 

If a^N were too large, DM particles in a galactic halo 
would scatter too frequently with the baryonic disk of a 
spiral galaxy, and would significantly disrupt it. Using 
the integrity of the Milky Way disk, Starkman et al. 
restrict the cross section to a^.^ < 5 x 10^^'*(m^/GeV) 
cm^. A more detailed study by Natarajan et al. f§\ re- 
quires a-^N < 5 X 10~'^^{m^/GeV) cm^. Both of these 
limits consider DM scattering only with hydrogen. As 
shown below in Eq. (|12p. the spin-independent DM- 
nucleon cross section scales as for large m^, and 
though the number density of helium {A = 4) is about 
10 times less than that of hydrogen {A — 1), taking it 
into account could improve these constraints by ~ 256/10 
~ 25. Chivukula et al. [40| showed that charged dark 
matter could be limited through its ionizing effects on 



interstellar clouds; this technique could be adapted for 
strongly interacting dark matter. 

Strong scattering of DM and baryons would also af- 
fect the cosmic microwave background radiation. Adding 
stronger DM-baryon interactions increases the viscos- 
ity of the baryon-photon fluid A strong cou- 
pling of baryons and DM would generate denser clumps 
of gravitationally-interacting matter, and the photons 
would not be able to push them as far apart. The peaks in 
the cosmic microwave background power spectrum would 
be damped, with the exception of the first one. The re- 
sulting constraint is a^N < 3 x 10~^'*(r7i^/GeV) cm^ 
and is not shown in Fig. [TJ These results do take helium 
into account, but do so only using instead of A'^. This 
possible change, along with the much more precise cosmic 
microwave background radiation data available currently, 
calls for a detailed re- analysis of this limit, which should 
strengthen it. 

Cosmic ray protons interact inelastically with interstel- 
lar protons, breaking the protons and creating neutral pi- 
ous that decay to high-energy gamma rays. A similar sit- 
uation could occur with a cosmic ray beam on DM targets 
instead @. The fundamental interaction is between the 
quarks in the nucleon and the DM; it is very unlikely that 
all quarks will be struck equally, and the subsequent de- 
struction of the nucleon creates pions. If the DM-nucleon 
cross section were high enough, the resulting gamma rays 
would be readily detectable. Prom this, Cyburt et al. Q 
place an upper limit of a^N < 7.6 x 10~^^(m^/GeV) 
cm^. Improvements could probably be made easily with 
a more realistic treatment of the gamma-ray data. 



B. Direct Detection Constraints 

Underground detector experiments have played a large 
role in limiting DM that can elastically scatter nuclei, 
giving the nuclei small but measurable kinetic energies. 
Due to the cosmic ray background, this type of detector is 
located underground. The usual weakly interacting DM 
candidates easily pass through the atmosphere and Earth 
en route to the detector. However, for large a^N the DM 
would lose energy through scattering before reaching the 
detector, decreasing detection rates. 

Albuquerque and Baudis have explored constraints 
at relatively large cross sections and large masses us- 
ing results from CDMS and EDELWEISS. In Fig.[Tl we 
present a crude estimate of the current underground de- 
tector exclusion region. The top line is defined by the 
ability of a DM particle to make it through the atmo- 
sphere [4l| and Earth to the detector without losing too 
much energy f7l| . The lower left corner and nearby points 
are taken from the official CDMS papers 6] with the 
aid of their website [l^]. The right edge is taken from 
DAM A 0. As the mass of the DM increases, the num- 
ber density (and hence the flux through Earth) decreases. 
At the largest values, the scattering rate within a fi- 
nite time vanishes. Finally, we have extrapolated each 
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of these constraints to meet each other, connecting them 
consistently. We call for a complete and official analysis 
of the exact region that CDMS and other direct detectors 
exclude. Our focus is on the cross sections in between the 
underground detectors and astrophysical limits. 

To investigate cross sections in this middle range, di- 
rect detectors must be situated above Earth's atmo- 
sphere, in high-altitude balloons, rockets, or satellites. 
Several such detectors have been analyzed for this pur- 
pose, though they were not all originally intended to 
study DM. Since these large (j^m limits have in some 
cases been calculated by people not connected with the 
original experiments, some caution is required. Never- 
theless, in Fig. [T] we show the claimed exclusion regions, 
following Starkman et al. [5| and Rich et al. along 
with Wandclt et al. [ll| and Erickcek et al. (including 
the primary references [H, |U EE 13 ) ■ We are primar- 
ily in accordance with Erickcek et al. These regions span 
masses of almost 0.1 GcV to 10^^ GeV, and cross sections 
between roughly 10^^'^ cm^ and 10"^^ cm^. These in- 
clude the Pioneer 11 spacecraft and Skylab, the IMP 7/8 
cosmic ray silicon detector satellite, the X-ray Quantum 
Calorimetry experiment (XQC), and the balloon-borne 
IMAX. These regions are likely ruled out, but not in ab- 
solute certainty, and there are gaps between them. The 
Pioneer 11 region is completely covered by the IMP 7/8 
and XQC regions, and is therefore not shown in Fig. [1] 
The region labeled RRS is Rich et al.'s analysis of a sil- 
icon semiconductor detector near the top of the atmo- 
sphere, truncated according to Starkman et al., and ad- 
justed with the appropriate ^-scaling as in Eq. 



III. EARTH'S HEAT FLOW 

Heat from the Sun warms Earth, but it is not retained. 
If all the incident sunlight were absorbed by Earth, the 
heating rate would be about 170,000 TW [2^. Some of it 
is reflected by the atmosphere, clouds, and surface, and 
the rest is absorbed at depths very close to the surface 
and then re-radiated Earth's blackbody tempera- 

ture would be about 280 K, and it is observed to be be- 
tween 250 and 300 K, supporting the idea of Earth-Sun 
heat equilibrium. Internal heating therefore has minimal 
effects on the overall heat of Earth 31 1 . 

Our focus is on this internal heat flow of Earth, as mea- 



sured underground. Geologists have extensively studied 
Earth's internal heat for decades [l^]. To make a mea- 
surement, a borehole is drilled kilometers deep into the 
ground. The temperature gradient in that borehole is 
recorded, and that quantity multiplied by the thermal 
conductivity of the relevant material yields a heat flux 

[13, El]. 

The deepest borehole is about 12 kilometers, which is 
still rather close to Earth's surface. Typical tempera- 
ture gradients are between 10 and 50 K/km, but these 
cannot hold for lower depths. If they did, all rock in 
the deeper parts of Earth would be molten, in contra- 



diction to seismic measurements, which show that shear 
waves can propagate through the mantle Current 
estimates place temperature gradients deep inside Earth 
between 0.6 and 0.8 K/km 

More than 20,000 borehole measurements have been 
made over Earth's surface. Averaging over the conti- 
nents and oceans, there is a heat flux of 0.087 ± 0.002 
W/m^ j23l Wf^. Integrating this flux over the surface of 
Earth gives a heat flow of 44.2 ± 1 TW [13, E^. Again, 
the heat flux is directly measured underground, all over 
Earth, and is independent of the solar flux. Earth's at- 
mosphere, and anything else above Earth's surface. Ob- 
viously, the possibility to make direct heat flow measure- 
ments under the surface is unique to Earth. 

While the heat flow value is known well, the origin 
of the heat is not, and in fact is undergoing much the- 
oretical debate [H, [sO]- Some specific contributors are 
known, however. The decay of radioactive elements pro- 
duces a significant amount; uranium and thorium de- 
cay in the crust generates about forty percent of the to- 
tal (Ssj. Potassium adds to this, though there is much 
less of it in the crust. However, there is potentially a 
large amount in the mantle and perhaps even the outer 
core [2^. KamLAND has a hint of detected neutri- 
nos coming from uranium and thorium decays (5l| . and 
it (along with other detectors) could potentially help 
to make the heat contribution from them more accu- 
rate [13, m 111 y, m, 113, S m, S mi . Larger concen- 
trations of uranium and thorium are excluded by Kam- 
LAND, and theoretical predictions from the Bulk Silicate 
Earth model are consistent with the forty percent value 
[55I IsgI . Wf\ ■ The remaining heat is due to processes in 
the core and perhaps even the mantle, although specific 
knowledge of Earth's interior is limited [B^I ■ 

The residual heat flow, which we assume to be 20 
TW [4^, [50| , we use as the target limit for the heat flow 
from DM annihilation. Models give values of the core's 
heat output between 2.3 TW and 21 TW, supporting the 
conservative choice of 20 TW Annihilation scenar- 
ios creating heat flows greater than 20 TW are therefore 
excluded. In fact, if heating by DM annihilation is im- 
portant at all, we show that it typically would exceed this 
value by more than two orders of magnitude. It is im- 
portant to note that we are not trying to solve geological 
heat problems with DM, and in fact our analysis implies 
it is very unlikely that DM is contributing to Earth's 
internal heat flow, which is interesting in itself. 



IV. DARK MATTER CAPTURE RATE OF 
EARTH 

The DM mass density, — n^m^^ in the neighbor- 
hood of the solar system is about 0.3 GeV/cm^ [l|. Nei- 
ther the mass nor the number density are separately 
known. The DM is believed to follow a nonrelativistic 
Maxwell-Boltzmann velocity distribution with an aver- 
age speed of about 270 km/s. If a DM particle scatters a 
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sufficient number of times while passing through Earth, 
its speed wih fall below the surface escape speed, 11.2 
km/s. Having therefore been gravitationally captured, it 
will then orbit the center of Earth, losing energy with 
each subsequent scattering until it settles into a ther- 
mal distribution in equilibrium with the nuclei in the 
core. For the usual weak cross sections Earth is effec- 
tively transparent, and scattering and capture are very 
inefficient. In contrast, we will consider only large cross 
sections for which capture is almost fully efficient. Note 
that for our purposes, the scattering history is irrelevant 
as long as capture occurs; in particular, the depth in the 
atmosphere or Earth of the first scattering has no bear- 
ing on the results. The energies of the individual struck 
nuclei are also irrelevant, unlike in direct detection ex- 
periments. We just require that the DM is captured and 
ultimately annihilated. 



from the atmosphere [23, l29| ; note also that the exclu- 
sion region in Section 4 would be unaffected by taking 
this into account, since the DM heating of Earth would 
still be excessive. 

The type of nucleus with which DM scatters depends 
on its initial trajectory through Earth. For a minimum 
path length of 0.2 R^, this trajectory runs through the 
crust, where the density is 3.6 g/cm'^ [63i] . and the most 
abundant element is oxygen [l6[. Choosing this path 
length and density are conservative steps. Any larger 
path length would result in more efficient capture, and 
a higher density and heavier composition (correspond- 
ing to a larger chord and therefore a different target nu- 
cleus, such as iron, which is the most abundant element 
in the core) would as well. A more complex crust or 
mantle composition, such as 30% oxygen, 15% silicon, 
14% magnesium and smaller contributions from other el- 
ements 



16j . would stop DM ^-^2 times more effectively. 



A. Maximum Capture Rate 



We begin by considering the maximum possible cap- 
ture rate of DM in Earth, which corresponds to Earth 
being totally opaque. Although our final calculations 
will involve the microscopic scattering cross section of 
DM on nuclei, this initial example deals with just the 
macroscopic geometric cross section of Earth. The flux 
per solid angle of DM near Earth is ny^v^/Air, where 
is the DM number density, and is the average DM ve- 
locity. Since Earth is taken to be opaque, the solid angle 
acceptance at each point on the surface is 27r sr. Thus 
the flux at Earth's surface is The capture rate 

is then found by multiplying by Earth's geometric cross 
section, = 47ri?^ ~ 5.1 x 10^^ cm^. Since is not 
known, this is {p-^/m^)aQV^. For — 270 km/s, this 
maximal capture rate is 



^niax ^ 2 X 10^^ 



GcV 



(1) 



We will show that our results depend only logarithmi- 
cally on the DM velocity, and hence arc insensitive to 
the details of the velocity distribution. 

This maximal capture rate estimate is too simplistic, as 
it assumes that merely coming into contact with Earth, 
interacting with any thickness, will result in DM cap- 
ture. Instead, we define opaqueness to be limited to path 
lengths greater than 0.2 R^, a value that incorporates the 
largest 90% of path lengths through Earth. This reduces 
the capture rate, but only by about 2%. We therefore 
adopt 0.2 i?0 as our minimum thickness to determine 
efficient scattering. This length, translated into a chord 
going through the spherical Earth, defines the new effec- 
tive area for Earth. The midpoint of the chord lies at 
a distance of 0.99Rq from Earth's center. Thus, prac- 
tically speaking, nearly all DM passing through Earth 
will encounter sufficient material. The above require- 
ments exclude glancing trajectories from consideration, 
for which there would be some probability of reflection 



B. Dark Matter Scattering on Nuclei 

When a DM particle (at ~ lO^^c) elastically scat- 
ters with a nucleus/nucleon (at rest) in Earth, it de- 
creases in energy and velocity. After one scattering with 
a nucleus of mass m^, DM with mass and initial 
velocity Vi will have a new velocity of 



Jl-2 , "^^"^"., -(l-cosg,„0, (2) 
1 _ 2^(1- cos 0,™). (3) 



All quantities are in the lab frame, except the recoil angle, 
^cm, which is most usefully defined in the center of mass 
frame (see Landau and Lifschitz 64]). Here and below 
we give the large limit for demonstration purposes, 
but we use the full forms of the equations for our results. 
After scattering, the DM has a new kinetic energy. 



KE- = 
1 



2 



1 



2!!^(l_cos0, 



The nucleus then obtains a kinetic energy of 



(5) 



KE^ = KE^^ -KE^ 



TO^TO^ 



(to^ -I- TO^)' 



:(1 



KE^'' 2^(l-COS0er„) 



TOx 



■mAv1{l - COB 9 cm) 



(6) 
(7) 
(8) 



6 



From the kinetic energy, the momentum transfer in the 
large hmit is: 



KE 



- niAvf {I- COS e cm) (9) 

(10) 



2mA 



In order to maintain consistency with others, we work 
with n and a in nucleon units even though the target we 
choose (oxygen) is a nucleus. This means that ua (where 
A represents the mass number of the target) is 



Ua 



A wlmA 



(11) 



In turn, the cross section for spin-independent s-wave 
elastic scattering is represented as 



XN 



(12) 



A'^a^N- 



Here A is the mass number of the target nucleus, which 
equals mA/mN, and fj,{A or N) is the reduced mass of 
the DM particle and the target. 

The A^ factor arises because at these low momentum 
transfers, the nucleus is not resolved and the DM is as- 
sumed to couple coherently to the net "charge" - the 
number of nucleons. (If this coherence is somehow lost, 
a factor A would still remain for incoherent scattering.) 
The momentum transfer q ~ y^2mAKE^ ~ ruAVi cor- 
responds to a length scale of ~ 10 fm for oxygen, much 
larger than the nucleus. We find that the corresponding 
nuclear form factor when the DM mass is comparable to 
the target mass is ~ 0.99. The square of the reduced mass 
ratio arises from the Born approximation for scattering, 
which is based on the two-particle Schrodinger equation 
cast as a singleparticle with relative coordinates and re- 
duced mass The spin-dependent scattering cross 
section does not have the A'^ factor in Eq. (fT2|) [3]. Our 
constraints could be scaled to represent this case by also 
taking into account the relative abundance of target nu- 
clei with nonzero spin in Earth, which is of order 1%. 

Note that if m^^ = ruA, and 9cm = t^, the DM can 
transfer all of its momentum to the struck nucleus, losing 
all of its energy in a single scattering through this scatter- 
ing resonance [16]. Taking this into account would make 
our constraints stronger over a small range of masses, but 
we neglect it. The nuclear recoil energy from this reso- 
nance is then ^m^vf. Since Vi is on average 270 km/s, 
this means that the maximum energy transferred from a 
collision is ^ 10"^ that of the annihilation energy, niyC^. 



C. Dark Matter Capture Efficiency 

From the full or approximate form of Eq. ([5]), we see 
that the DM kinetic energy is decreased by a multiplica- 



tive factor that is linear in cos 0cm- If, in each indepen- 
dent scattering, we average over cos 0cm, the average fac- 
tor by which the kinetic energy is reduced in one or many 
scatterings will simply be that obtained by setting cos 0cm 
— throughout. (For s-wave scattering, the cos 0cm dis- 
tribution is uniform.) 

We will define efficient capture so that the heating is 
maximized. To be gravitationally trapped, a DM particle 
must be below the escape speed of Earth (t^csc = 11-2 
km/s), or equivalently, its kinetic energy must be less 



than ■ 



After one scattering event, the DM kinetic 



energy is reduced: 

KE/'^KE^'^fim^). (13) 
In successive collisions, this is compounded until 



(14) 



Note that for coUinear scatterings, the velocity loss in 
Eq, ^ is also speed loss, leading to the same definition 

of iVscat- 

Therefore, on average, the number of scatterings re- 
quired to gravitationally capture the DM is 



-2 ln(wi/i;csc) 



In 



mAm^ 



— In (Vi/Vcsc), 
TTIA 



(15) 
(16) 



where we have set cos 0cm = 0, since this corresponds 
to the average fractional change in the kinetic energy. 
Again, for simplicity the same element is taken to be the 
target each time. Note that since the initial DM velocity 
is inside the logarithm, A^scat is insensitive to even large 
changes in the assumed initial velocity. 

The number of scatterings for a given mass is large. A 
DM particle that has the same mass as the target nucleus 
will scatter about 10 times before it is captured. Note 
that the required A'scat scales as in the large mass 
limit, becoming very large: for above 16 TeV (10'^ 
times the target mass), A'scat is already larger than 3000. 
The actual energy losses in individual collisions are irrel- 
evant for our analysis, as we require only that the DM is 
captured after many collisions. For large values of A'scat, 
all scattering histories will be well-characterized by the 
average case. 

So far, these equations have just been kinematics; the 
required A'scat for stopping has not yet been made specific 
to Earth. It becomes Earth-specific by relating A'scat to 
the path length in Earth L and the mean free path A: 



L 

Nscat = ^ = LUACT^A- 



(17) 



The column density of Earth then defines the required 
cross section to generate A'scat scatterings. The short- 
est path the particle could travel in is a straight line, so 
we use that as the minimum. Any other path would be 
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longer, and hence more effective at capture. Tliis there- 
fore defines the most conservative limit on ct^a- Since we 
have fixed cos ^cm to be on average, in fact the path 
will not be completely straight. However, the lab frame 
scattering angles are small. 

For elastic collisions between two particles, the range 
of scattering angles in the lab frame depends on the two 
masses, mi and m2. There is a maximum scattering an- 
gle when one mass is initially at rest in the lab frame (in 
this case, 7712) [64|. Ifwi < 7712, there is no restriction on 
the scattering angle, which is defined in relation to mi's 
initial direction (m2 is at rest). However, if mi > m2, 
then 



Our main focus is 



sm( 



mi 



lab 



m2/mi. 



(18) 



> m2 = iriA- For m^ some- 



what greater than m^, note that the DM scattering angle 
in the lab frame is always very forward. 

Combining Eqns. ([12]), ([IS]), and ([TT]) . the minimum 
required cross section to capture a DM particle is 



rriA 



pL 



-^scat("7x), 



(19) 



-2 ln(i;i/wosc) 



In 



\ — 2 ""^""x 

(nix+mA)'^ 



\mA J pL 



Again, we choose a path length of 0.2 i?^, to select about 
90% of the path lengths in Earth. Taking this length as 
a chord through Earth, the location corresponds to the 
crust, with an average density of 3.6 g/cm"^, where the 
most common element is oxygen. We also choose an in- 
coming DM velocity of 500 km/s, which effectively selects 
the entire thermal distribution. A slower DM particle is 
more easily captured. These parameters give a required 
cross section of 



'1.8 X 10" 



^(1) 

m(16) 



In 



-| rj 16 GcV 7yix 
^ (m^ + ie GeV)^ 



2.2 X 10"^"^ cm' 



(20) 
(21) 



Note that we use the unapproximated version, Eq. ([^0]) . 
for our figure, and give the large limit in the equations 
for demonstrative purposes. When is comparable to 
m^, cr™j^J^ is different from the approximated, large m^ 
case in an important way. 

The resulting curve for cr™^^ is shown in Fig. |2l as 
the lower boundary of the heavily-shaded exclusion re- 
gion. The straight section of this constraint is easily seen 
from Eq. ([2T|) . as the required cross section for our ef- 
ficient capture scenario scales as m^, due to the large 
number of collisions required for stopping, as in Eq. ([H 




FIG. 2: Inside the heavily-shaded region, dark matter anni- 
hilations would overheat Earth. Below the top edge of this 
region, dark matter can drift to Earth's core in a satisfac- 
tory time. Above the bottom edge, the capture rate in Earth 
is nearly fully efficient, leading to a heating rate of 3260 TW 
(above the dashed line, capture is only efficient enough to lead 
to a heating rate of > 20 TW). The mass ranges are described 
in the text, and the light-shaded regions are as in Fig. [1] 



At lower masses, the curved portion has its minimum 
at the mass of the target. DM masses close to that of 
the target can be captured with smaller cross sections 
because a greater kinetic energy transfer can occur for 
each collision. At very low masses, much less than the 
mass of the target, the DM mass dependence in the log- 
arithm is approximated differently. In this limit, a^^^jsi is 
~ 10"-^^ (GeV/m^). As the DM mass decreases, it be- 
comes increasingly more difhcult for the DM to lose en- 
ergy when it strikes a nucleus. As noted above, the cross 
section constraints in the spin-dependent case could be 
developed, and would shift the results up by 3 or 4 orders 
of magnitude. All of the other limits that depend on this 
coherence factor would also shift accordingly. 
The lower edge of the exclusion region is generally 
rather sharp, because of these parameters. For example, 
consider the case of large m^, where A^scat is also large. If 
the corresponding cross section is decreased by a factor 5, 
so is the number of scatterings, and by Eq. (fT4]) . the com- 
pounded fractional kinetic energy loss would only be the 
1/(5 root of that required for capture. For small cross sec- 
tions, as usually considered, the capture efhciency is very 
low. To efficiently produce heat, the minimum cross sec- 
tion must result in ~ 90% DM capture. We stress again 
that we are not concerned with where the DM is captured 
in Earth, so long as it is. The probability for capture can, 
however, be decreased using Poisson statistics (shown in 
Fig. [2] as the dashed line with the accentuated dip at low 
masses) to yield just 20 TW of heat flow. This extension 
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and the upper edge of the exclusion region are described 
below. 



V. DM ANNIHILATION AND HEATING 
RATES IN EARTH 

A. Maximal Annihilation and Heating Rates 

Once it is gravitationally captured, DM will continue 
to scatter with nuclei in Earth, losing energy until drift- 
ing to the core. Once there, because of the large cross 
section, the DM will thermalize with the nuclei in the 
core. The number of DM particles 9i is governed by the 
relation between the capture {Tc) a-nd annihilation (F^,) 
rates 1661: 



1 



1 



(22) 



We neglect the possibility of evaporation [15( for the mo- 
ment, which will affect our results for low and ctj^at, 
as we will explain further below. The variable t is the age 
of the system, is related to the DM self-annihilation 
cross section tr^^ by 

(23) 



off 



where VcS is the effective volume of the system '66'] . For 
the relevant cross sections considered, equilibrium be- 
tween capture and annihilation is generally reached (see 
below), so the annihilation rate is 



F^ = -Fc 
2 



(24) 



The effective volume is determined by the method of 
Griest and Seckel (1988) [sB], which is essentially the vol- 
ume of the DM distribution in the core. The number den- 
sity of DM is assumed to be an exponentially decaying 
function, exp(— m^(/)/fcr), like the Boltzmann distribu- 
tion of molecules in the atmosphere. The temperature of 
the DM in thermal equilibrium is T . The variable <j) is 
the gravitational potential, integrated out to a radius r, 
written as 

GM{f) 



-df\ 



M{f) 



(25) 



(26) 



47r / r'^p{r')dr' . 
Jo 

The resulting effective volume using the radius of Earth's 
outer core, approximately O ARff, , a temperature of 5000 
K, and a density of 9 g/cm'^ ^3], is 



Vm 



An 



r e dr. 



(27) 



1.2 X 10^^ cm^ 



100 GcV 



0.5\ 



i^e"" du 







24 o/100GeV\^ 
5.3 X 10^^ cm^ 



(28) 



For increasing , the integral (without the prefactor) in 
the second line of Eq. ([?7|) quickly reaches an asymptotic 
value of about 0.44. 

At very large masses, the effective volume for annihi- 
lation becomes very small. For instance, at > 10^° 
GeV, the radius of the effective volume is < 0.1 km. With 
such a large rate of energy injected in such a small vol- 
ume, the core temperature would be increased, requiring 
a more careful treatment. However, in Section 5B, we 
state how a limit on the annihilation cross section from 
the unitarity condition 67, 68, 69] makes us truncate our 
bound at my. = 10^" GeV, as reflected in Fig. [21 There- 
fore, this small effective volume is not a large concern for 
our exclusion region. 

We assume that the DM annihilates into primarily 
Standard Model particles, which will deposit nearly all 
of their energy into Earth's core (with small corrections 
due to particle rest masses and the escape of low-energy 
neutrinos). When all of the DM captured is efficiently 
annihilated, as specified, the heating rate of Earth is in 
equilibrium with the capture rate: 

Fhcat = Fc X = a^s 



Px 



27r(0.99i?e)^(270 km/s) (29) 



= 3260 TW. 

This heat flow is independent of DM mass, since the flux 
(and capture rate, when capture is efficient) scales as 
l/niy, while each DM particle gives up in heat when 
it annihilates. The value is much larger than the mea- 
sured rate of 44 TW we discussed in Section 3. 



B. Equilibrium Requirements 

Does the timescale of Earth allow for equilibrium be- 
tween capture and annihilation for our scenario? In order 
for Eq. to be in the equilibrium limit, taiih^ (ty/Tc^) 
must be of order unity. This is true if t\/Tc^ has a value 
of a few or greater. Since Earth is about 4.5 Gyr old, 
we conservatively require that the time taken to reach 
equilibrium should be less than about 1 Gyr. From this, 
a realistic annihilation cross section is found. The condi- 
tion 



tanh^(tv/f^) = 1; ty/Tc^ =i few (30) 
allows the relation 



> 



(few)2 



(31) 
(32) 



For an efficient capture rate (Eq. [T|), the time of 1 Gyr, 
and the limit of large m^, this requires an annihilation 
cross section for equilibrium of 



/r -\r\ 



V ^x 



cmVs. (33) 
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Since this required lower bound is much smaller than 
that of typical weakly interacting DM particles that are 
thermal relics (< cr^^v >~ 10^^^ cm^/s [l|), it should be 
easily met. One expects large scattering cross sections to 
be accompanied by large annihilation cross sections, so 
that even the possibility of p-wave-only suppression of 
the annihilation rate should not be a problem. 

For very large masses, the required annihilation cross 
sectionm, while small, approaches a quantum mechanical 
limit. For example, for > 10^° GeV the s-wave cross 
section exceeds the unitarity bound ^67|, (68,, £9J • We note 
that this may also affect constraints on supermassive DM 
based on neutrinos from annihilations [sj, [H, [13] . To be 
conservative, wc therefore do not extend our constraints 
beyond this point, though they may still be valid. 

The timescale also has to be long enough for DM to 
drift down to the core. If a-)^N is too large, the DM 
will experience too many scatterings and will not settle 
into the core, and thus may not annihilate efficiently. 
Following Starkman et al. we define the upper edge 
of our exclusion region to require a drift time of < 1 Gyr. 
This places a restriction of 



< 7.7 X lO^^^cm^ 



(my/GeV) 
A^p{A)/^i{N)y 



< 7.7 X 10-^"cm 



20„_2 mJGeV 



^4 / mx+™w y 



a^M < 2.5 X 10-23 cm^ (^) , (34) 



for a target of iron. However, a more detailed calcula- 
tion might relax this requirement. For example. Stark- 
man et al. j2.] show that for large values of the capture 
cross section and certain other conditions, annihilation 
may be efficient enough to occur in a shell, before the 
DM reaches the core. This would generally still be sub- 
ject to our constraint on heat from DM annihilation, and 
hence our exclusion region might extend to larger cross 
sections than shown. The two features of this drift line 
at low mass occur around the mass of the target and the 
mass of a nucleon, due to the various dominances of the 
mass-dependent term in the denominator. The details of 
the shape of this drift line at low masses are irrelevant, 
because the astrophysical constraints already exclude the 
corresponding regions. 

Aside from drifting to the core, the question of whether 
heavy DM can actually get to Earth has been asked [tO, 
[Tlj . The low- velocity tail of the high-mass DM thermal 
distribution in the Solar System may be driven into the 
Sun by gravitational capture processes |70i,i71|, especially 
because this DM's velocity is on the order of the orbital 
speed of Earth in the Solar System, which is about 30 
km/s. However, this would affect only a tiny fraction of 
the full thermal distribution that we require to be effi- 
ciently captured. 



C. Annihilation and Heating Efficiencies 

We are not picking a specific model for the annihilation 
products, aside from considering only Standard Model 
particles, which will deposit their energy in Earth, with 
the exception of neutrinos. Our constraint thus has a 
very broad applicability. As noted the calculated heat 
fiow if DM annihilation is important is 3260 TW, which 
is very large compared to our adopted limit on an uncon- 
ventional source of 20 TW (or even the whole measured 
rate of 44 TW). Typically then, either DM annihilation 
heating is overwhelming or it is negligible, inside or out- 
side of the excluded region. As shown in Section 4, the 
kinetic energy transferred from DM scattering on nuclei 
is about 6 orders of magnitude less than the energy from 
DM annihilations. This contribution to Earth's heat is 
too low to be relevant for global considerations. However, 
it would be interesting to consider the more localized ef- 
fect of the kinetic energy deposition in the atmosphere 
for very large cross sections. 

There are circumstances in which the heating from DM 
annihilations can take a more intermediate value, includ- 
ing down to the chosen 20 TW number. As explained 
above, typically the number of scatterings required to 
gravitationally capture the DM is very large. Therefore, 
a small decrease in (T^^,^r and the proportionate change in 
the expected number of scatterings means that the com- 
pounded kinetic energy loss is nearly always insufficient. 
However, at low m^, the number is small enough that 
upward fluctuations relative to the expected number can 
lead to capture. If Ngcat collisions typically lead to effi- 
cient capture for a cross section cr™^\ as defined above, a 
new and smaller N may be defined by the condition that 
the Poisson probability Prob(A^ > iVscat) = 20 TW / 
3260 TW = 1/163. With this TV, and its proportionately 
smaller a^N , upward Poisson fluctuations in the number 
of scatterings lead to efficient capture for a fraction 1/163 
of the incoming flux. Note that this small capture frac- 
tion is not just the low-velocity tail of the DM thermal 
distribution, since we have deflned these conditions for 
the highest incoming velocities, Vi = 500 km/s. 

The resulting constraint on cr^jv is shown by the dashed 
line that dips below the main excluded region in Fig. [21 
The enhanced valley around 16 GeV again arises from 
the ease of capture when the DM mass is near the target 
mass. Note that for each mass the required < Cxx*^ ^ 
increased by the same factor that decreased the original 
required cr™^. Since most of this exclusion region is al- 
ready covered by underground detectors, its details may 
not be so important. 

For low DM masses, evaporative losses of DM from the 
core due to upscattering by energetic iron nuclei may be 
relevant [lB| • Simple kinematic estimates show that DM 
masses below ~ 5 GeV might be affected. However, this 
is only potentially important if a-^^N is small enough that 
scattering is very rare - since otherwise any upscattered 
DM will immediately downscatter. From the considera- 
tions above about Poisson fluctuations in the number of 
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scatterings, we expect that this should only be relevant 
between the dark-shaded region and the dashed line. 



VI. DISCUSSION AND CONCLUSIONS 
A. Principal Results 

As summarized in Fig. [1] while very large DM-nucleon 
scattering cross sections are excluded by astrophysical 
considerations, and small cross sections are excluded by 
underground direct DM detection experiments, there is 
a substantial window in between that has proven very 



difhcult to test, despite much effort [ 2 | .ll0l. llll . ll2l ll: 

11, [13, E S Us Hi, M m, M MmMmm 

735. High-altitude experiments have excluded only parts 
of this window. In this window, DM will be efRciently 
captured by Earth. We point out that the subsequent 
self-annihilations of DM in Earth's core would lead to 
an enormous heating rate of 3260 TW, compared to the 
geologically measured value of 44 TW. 

We show that the conditions for efRcient capture, anni- 
hilation, and heating are all quite generally met, leading 
to an exclusion of a^N over about ten orders of mag- 
nitude, which closes the window on strongly interact- 
ing DM between the astrophysical and direct detection 
constraints. These new constraints apply over a very 
large mass range, as shown in Fig. [51 We have been 
quite conservative, and so very likely an even larger re- 
gion is excluded. These results establish that DM in- 
teractions with nucleons are bounded from above by the 
underground experiments, and therefore that these in- 
teractions must be truly weak, as commonly assumed. 
This means that direct detection experiments are look- 
ing in the correct a-^N range when sited underground and 
motivates further theoretical study of weakly interacting 
DM [li, [zi, [73, [Zi- Furthermore, it means that DM- 
nucleon scattering cannot have any measurable effects in 
astrophysics and cosmology, and this has many implica- 
tions for models with strongly or moderately interacting 
DM 0, [zl, [tI, [tI] and other astrophysical constraints on 
the DM-nucleon interaction cross section J9g ^8Ql . This 
exclusion region also completely covers the cross section 
range in which strongly interacting dark matter might 
bind to nuclei [sil]. 

To evade our constraints, extreme assumptions would 
be required: that DM is not its own antiparticle, or that 
there is a large ( > 163) particle-antiparticle asymme- 
try from DM, or that DM self-annihilations proceed only 
to purely sterile non-Standard Model particles, at the 
level of > 100:1. (Although in Ref. [631 such a large 
branching ratio to neutrinos was considered, it was em- 
phasized that this was used only to set the most conser- 
vative bound on the DM annihilation cross section, and 
not to be representative of a realistic model.) While our 
constraints are based on Earth's measured heat flow, it 
is important to emphasize that DM capture and anni- 
hilation generally cannot contribute measurably without 



being overwhelming, and hence are excluded. Thus, in 
the ongoing debate over the unknown sources of Earth's 
heat flow, it seems that DM can play no role. The 
most important next step in refining our understanding 
of the known generators of the measured 44 TW will 
come from isolating the contribution from uranium and 
thorium decays by measuring the corresponding neutrino 
fluxes i5^, iM, i55,, ^ i57„ iM, i59,, M, SL] ■ 



B. Comparison to Other Planets 

Starkman et al. [2] calculated the efficient Earth cap- 
ture line for DM, but only had model-dependent results. 
We have now considered the consequences of annihila- 
tion in Earth, and have shown that it gives a model- 
independent constraint. Other planets have been dis- 
cussed, such as Jupiter and Uranus [H, HI, US IH, 
i3, m, m, [131, but Earth is the best laboratory. It 
is the best understood planet, with internal heat flow 
data measured directly underground, from many loca- 
tions, and Earth's composition and density profile are 
weU known [2l[43,|48, 49, 50, 62, 63]. Importantly, the 
relative excess heat due to DM annihilation would be 
much greater for Earth than the Jovian planets. 

What about other planets? The maximal heating rate 
due to DM scales with surface area, and can be compared 
with the internal heating rates estimated from infrared 
data [3l| . If a constraint can be set, the minimum cross 
section cr™j\f that can be probed scales with the planet's 
column density as {nL)~^ (see Eq. (ITTl) ') up to nontrivial 
corrections for composition (see Eq. (fT2)) ). Note that the 
column density nL is proportional to the surface grav- 
ity ^ GM/E?^ which varies little between the planets, 
as noted in Table II [31]. Due to its known (and heavy) 
composition and well-measured (and low) internal heat, 
the strongest and most reliable constraints will be ob- 
tained considering Earth. As an interesting aside, it may 
then be unlikely that heating by DM could play a signifi- 
cant role in explaining the apparent overheatin g of some 
extra-solar planets ( "hot Jupiters" ) [H, [H, [s^ l85l. [s^ . 



C. Future Directions 

While it appears that the window on strongly inter- 
acting DM is now closed over a huge mass range, more 
detailed analyses are needed in order to be absolutely 
certain. Our calculations are conservative, and the true 
excluded region is likely to be larger. It would be espe- 
cially valuable to have new analyses of the astrophysical 
limits and the underground detector constraints. This 
would give greater certainty that no sliver of the window 
is still open. For DM masses in the range 1 GeV to 10^° 
GeV, the upper limits on the DM scattering cross section 
with nucleons from CDMS and other underground exper- 
iments have been shown to be true upper limits. Thus 
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TABLE 11: Comparison of potential dark matter constraints 
using various planets [sj. A greater difference between dark 
matter and internal heating rates give greater certainty. The 
minimum cross section probed scales roughly with the sur- 
face gravity. Earth is the best for setting reliable and strong 
constraints. 



Planet 


DM Max. Heating 


Internal Heat 


Surface Gravity 




(TW) 


(TW) 


(units of ®) 


Earth 


3.3 X 10^ 


0.044 X 10^ 


1.00 


Jupiter 


420 X 10^ 


400 X 10^ 


2.74 


Saturn 


290 X 10^ 


200 X 10^ 


1.17 


Uranus 


53 X 10^ 


< 10'^ 


0.94 


Neptune 


50 X 10^ 


3 X 10^ 


1.15 



the DM does indeed appear to be very weakly interact- 
ing, and it will be challenging to detect it. 
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